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We consider the Affleclc-Dine baryogenesis comprehensively in the minimal supersymmetric standard 
model with gauge-mediated supersymmetry breaking. Considering the high temperature effects, we 
see that the Affleck-Dine field is naturally deformed into the form of the Q ball. In the natural 
scenario where the initial amplitude of the field and the A-terms are both determined by the non- 
renormalizable superpotential, we obtain only very a narrow allowed region in the parameter space 
in order to explain the baryon number of the universe for the case that the Q-ball formation occurs 
just after baryon number production. Moreover, most of the parameter sets suited have already 
been excluded by current experiments. We also find new situations in which the Q-ball formation 
takes place rather late compared with baryon number creation. This situation is more preferable, 
since it allows a wider parameter region for naturally consistent scenarios, although it is still difficult 
to realize in the actual cosmological scenario. 
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I. INTRODUCTION 

The AfHeck-Dine (AD) mechanism [Q] is the most 
promising scenario for explaining the baryon number of 
the universe. It is based on the dynamics of a (complex) 
scalar field (j) carrying baryon number, which is called the 
AD field. During inflation, the expectation value of the 
AD field develops at a very large value. After inflation 
the inflaton field oscillates about the minimum of the ef- 
fective potential, and dominates the energy density of the 
universe like matter, while the AD field stays at the large 
field value. It starts the oscillation, or more precisely, ro- 
tation in its effective potential when H m^^eff, where 
H and m^^eff = are the Hubble parameter and 

the curvature of the potential of the AD field. Once it 
rotates, the baryon number will be created as ~ (^0^, 
where u) = 9 is the velocity of the phase of the AD field. 
When t ~ r^, where is the decay rate of the AD field, 
it decays into ordinary particles carrying baryon num- 
ber such as quarks, and the baryogenesis in the universe 
completes. 

However, important effects on the field dynamics were 
overlooked. It was recently revealed that the AD field 
feels spatial instabilities . Those instabilities grow very 
large and the AD field deforms into clumpy objects: Q 
balls. A Q ball is a kind of nontopological soliton, whose 
stability is guaranteed by the existence of some charge 
Q In the previous work we found that all the 
charges which are carried by the AD field are absorbed 
into formed Q balls, and this implies that the baryon 
number of the universe cannot be explained by the relic 
AD field remaining outside Q balls after their formation. 

In the radiation dominated universe, charges are evap- 
orated from Q balls , and they will explain the baryon 
number of the universe |^ . This is because the minimum 
of the (free) energy is achieved when the AD particles are 



freely in the thermal plasma at finite temperature. (Of 
course, the mixture of the Q-ball configuration and free 
particles is the minimum of the free energy at finite tem- 
perature, when the chemical potential of the Q ball and 
the plasma are equal to be in the chemical equilibrium. 
This situation can be achieved for very large charge of Q 
balls with more than 10^" or so Q.) Even if the radiation 
component is not dominant energy in the universe, such 
as that during the inflaton-oscillation dominant stage just 
after the infiation, high temperature effects on the dy- 
namics of the AD field and/or Q-ball evaporation are 
important. Therefore, in this article, we investigate the 
whole scenario of the AD mechanism for baryogenesis 
in the minimal supersymmetric standard model (MSSM) 
with the gauge-mediated supersymmetry (SUSY) break- 
ing in the high temperature universe. 

In Sec. II, we identify one of the flat directions as the 
AD field, and look for the form of the effective poten- 
tial in the context. We will note how the AD field pro- 
duces the baryon number in the universe in Sec. HI. In 
Sec. IV, we consider the dynamics of the (linearized) fluc- 
tuations. The Q-ball formations are simulated in Sec.V, 
and we will construct the formula of the Q-ball charge in 
terms of the initial amplitude of the AD held. Section 
VI has to do with the mechanism for the evaporation of 
the Q-ball charge, and we estimate the total charge evap- 
orated, which will be the baryons in the universe later. 
In Sec. VII, we seek for the consistent cosmological Q- 
ball scenario, and flnd it very difficult as contrary to our 
expectation. In Sec. VIII, we will investigate the Q-ball 
formation and natural scenario in more generic gauge- 
mediation model, where we take the different scales for 
the potential height and the messenger scale of gauge- 
mediation. In Sec. IX, we will consider a new situation 
where the baryon number creation takes place when the 
field starts the rotation, but the Q-ball formation occurs 
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rather later. In that case, we will find the wider consis- 
tent regions of the parameter space in some situations. 
The detections of the dark matter Q ball and their con- 
straints are described in Sec. IX. Section X is devoted for 
our conclusions. 



II. FLAT DIRECTIONS AS THE AFFLECK-DINE 
FIELD 

A flat direction is the direction in which the effective 
potential vanishes. There are many flat directions in 
the minimal supersymmetric standard model (MSSM), 
and they are listed in Rcfs. Since they consist of 

squarks and/or sleptons, they carry baryon and/or lep- 
ton numbers, and we can identify them as the Affleck- 
Dine (AD) field. Although the flat directions are exactly 
flat when supersymmetry (SUSY) is unbroken, it will be 
lifted by SUSY breaking effects. In the gauge-mediated 
SUSY breaking scenario, SUSY breaking effects appear 
at low energy scales, so the shape of the effective poten- 
tial for the flat direction has curvature of order of the 
electroweak mass at low scales, and almost flat at larger 
scales (actually, it grows logarithmically as the field be- 
comes larger). To be concrete, we adopt the following 
special form to represent such a kind of the potential : 



V{<P)=ml\og 




(1) 



where $ is the complex scalar field representing the flat 
direction. In spite of the specific form, it includes all 
the important features for the formation of Q balls, since 
they are formed at large field value, where the logarithmic 
function has to do with the dynamics of the field. [] 

Although we are considering the gauge mediation 
model, the flat directions are also lifted by the gravity 
mediation mechanism, since the gravity effects always ex- 
ist. We can write in the form 
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where the K term is the one-loop corrections pO| , and 
M = 2.4 X 10^^ GeV is the reduced Planck mass. Since 
the gravitino mass is much smaller than 1 TeV, 



^ In Ref. ^1, the effective potential has the form 

V ~ F^[log(02)]^ 

where F^^'^ >> m^. However, the dynamics is very similar 
to that in the potential and so is the Q-ball formation. 
When we seek for the allowed region in the parameter space, 
the difference between F^^^ and leads to the different 



which is suggested by the gravity-mediated SUSY break- 
ing scenario, this term will dominate over the gauge- 
mediation effects only at very large scales [||J^] . 

The flat directions may also be lifted by nonrenor- 
malizablc terms. When the superpotential is = 
A<I>"/(nM"-3), the terms 
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h.c. 



are added to the effective potential for the flat directions, 
where Va denotes A-terms and we assume vanishing cos- 
mological constant to obtain them, so \A\\ ~ ^3/2- 

In addition to the terms above, there are those terms 
which depends on the Hubble parameter _ff , during infla- 
tion and inflaton oscillation stage which starts just after 
inflation. These read as 



Vah — ^clhH— — 5- -I- h.c, 



(4) 
(5) 



where ch is a positive constant and an is a complex 
constant with a different phase from A\ in order for the 
AD mechanism to work; we need a large initial amplitude 
of the field and a kicking force for the phase rotation, 
which leads to the baryon number creation. 

Now we are going to consider thermal effects on the 
effective potential for the AD field. The fiat directions 
can be lifted by thermal effects of the SUSY breaking 
because of the difference of the statistics between bosons 
and fermions. If the AD field couples directly to parti- 
cles in thermal bath, its potential receives thermal mass 
corrections 
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where ct is proportional to the degrees of freedom of 
the thermal particles coupling directly to the AD field. 
Therefore, if the amplitude of the AD field is large 
enough, this term will vanish, because the particles cou- 
pled to the AD field get very large effective mass, nief / ^ 
/l^l :s> T, where / is a gauge or Yukawa couphng con- 
stant between those particles and the AD field, so they 
decouple from thermal bath. If we consider the forma- 
tion of large Q balls, the value of the AD field should 
be very large, and we can neglect this term in the most 
situations. 

On the other hand, there is another thermal effect on 
the potential. This effect comes from indirect couplings 
of thermal particles to the AD field, mediated by heavy 
particles which directly couple to the AD field, and it 
reads as 
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conclusion. See Sect. VIII 
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where 
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a^T) where a{T) = g\T)/ATT, 

estimated at the temperature T. We wiU combine these 
two terms in one form as 



Vt = r log 1 
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Therefore, comparing with Eq.(|i|), we can regard the ef- 
fective mass as 



m(T) 



7710 (T < m^) 
T (T > m^) 



and the effective potential as 



= m''(T) log 1 
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As can be seen later, the Q-ball formation takes place 
very rapidly, so the time dependence on T is less effec- 
tive, since T cx t^^/"* during inflaton-oscillation domi- 
nated universe. Therefore, we can regard T as constant 
and apply the numerical results obtained for ni{T) — 
to the case for m{T) = T hy rescaling the variables to be 
dimensionless with respect to T. 



III. AFFLECK-DINE MECHANISM 

It was believed that the Affleck-Dine (AD) mechanism 
works as follows During inflation, the AD fleld are 

trapped at the value determined by the following condi- 
tions: Vlf{<P) - V^jj{^) and VXh(^) ~ 0. Therefore, we 
get 



0™„^(i/Af"-3)i/("-2), 
sin[n9min + arg{aH)] ^ 0, 



(11) 



where ch, A ^ 1 are assumed, and $ ~ (0e'^)/-\/2. After 
inflation, the inflaton oscillates around one of the mini- 
mum of its cfl'cctive potential, and the energy of its os- 
cillation dominates the universe. During that time, the 
minimum of the potential of the AD field are adiabati- 
cally changing until it rolls down rapidly when ~ w = 
\V'{(j))/(l)\^/'^ . Then the AD field rotates in its potential, 
producing the baryon number ub — Ocj)^. After the AD 
field decays into quarks, baryogenesis completes. 

In order to estimate the amount of the produced 
baryon number, let us assume that the phases of A\ and 
an differ of order unity. Then the baryon number can be 
estimated as 
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where e = {m^j2/oj), n^'"^' = ujq 
in the second line. Notice that the contribution from 



and H 



is used 



the Hubble A-term is at most comparable to this. When 
e = 1, the trace of the motion of the AD field in the 
potential is circular orbit. If e becomes smaller, the orbit 
becomes elliptic, and finally the field is just oscillating 
along radial direction when e = 0. We call e as the 
ellipticity parameter below. 

When the logarithmic potential ( p^ ) is dominant, lo 
rn? /4>, so the ellipticity parameter is 



^3/200 

^'^ m2(T) ■ 
For £ ~ 1, it is necessary to have 

m(T) - (m;^/-iAf"-3)i/(2n-4)^ 



(13) 
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where we use Eq.(|TT|) with H ^ u ^ m'^{T)/(j)o. When 
m{T) ~ T, the conditions for the AD field to rotate cir- 
cularly are 
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GeV (n==6). 

(15) 



On the other hand, when the potential is dominated 
by Vgrav, ^ — "T-3/2- Then £ ~ 1 always holds. Of 
course, this conclusion does not necessarily true, if the 
A-terms do not come from the nonrenormalization term 
in the superpotential, such as in Eq.(H). e is generally 
arbitrary, if the origin of the A-terms is in the Kahler 
potential, which we have less knowledge. 



IV. INSTABILITIES OF AFFLECK-DINE FIELD 

In the Affieck-Dine mechanism, it is usually assumed 
that the field rotates in its potential homogeneously to 
produce the baryon number. However, the field feels spa- 
tial instabilities, and they grow to nonlinear to become 
clumpy objects, Q balls In this section, we study the 
instability of the field dynamics both analytically and 
numerically. 

For the discussion to be simple, we consider the field 
dynamics in the logarithmic potential only when the 
gauge-mediated SUSY breaking effects dominate. This 
is enough to investigate the whole dynamics, since fluc- 
tuations are developed when the field stays in this region 
of the effective potential. We write a complex field as 
$ — (0e*^)/'\/2, and decompose into homogeneous part 
and fiuctuations: (j) ^ (j) + Scj) and 9 ^ 6 + 66. Then the 
equations of motion of the AD field read as P,p^ 



3iJ0 - 9^ 



m?(f> 



= 0, 
= 0, 
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for the homogeneous mode, and 



- + 2 



- -^S<l) + V"{<l))S(t>^0, 

- 2^e6(f> - (l)^S9 = 0, (18) 



for the fluctuations, where a is the scale factor of the 
universe, and 
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where we assume </> ^ m in the last step. Equation (|lj 
represents the conservation of the charge (or number) 
within the physical volume: 9(^)^0? — const. 

Now let us forget about the cosmic expansion, since 
it only makes our analysis be complicated. Thus, a = 
const, here. It is natural to take a circular orbit of the 
field motion in the potential. (For a noncircular orbit, it 
is difhcult to discuss analytically, so we will show numer- 
ical results later.) Then we can take ip ^ (f>o = const., 
and the phase velocity 9 = {V ~ \/2mV0o- We 
seek for the solutions in the form 



5(f) = 6(j)Qe 
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(20) 



If a is real and positive, these fluctuations grow exponen- 
tially, and go nonlinear to form Q balls. Inserting these 
forms into Eqs.(18), we get the following condition for 
nontrivial (500 and 59q, 
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This equation can be simplified to be 
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= 0. (22) 



In order for a to be real and positive, we must have 
4m4 



a2 



fc2 

^ < 0. 



(23) 



Therefore, we obtain the instability band for the fluctu- 
ations as 



k 2m2 

a 00 



(24) 



We can easily derive that the most amplified mode ap- 
pears at {k/a)max = (3/2)^/2m2/0o, and the largest 
growth factor is a^ax = a.{kmax) = to2/(\/20o)- 

It is easier to decompose a complex field into its real 
and imaginary parts for analyzing the dynamics of the 



field when its motion is noncircular. For the numerical 
calculation, it is convenient to take all the variables to be 
dimensionless, so we normalize 'as (p — 0/to, k = k/m, 
T = mt, and ^ = mx. Writing ip — {(pi + i(^2)/\/2, we 
get the equations for the homogeneous mode as 
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^ ^0 
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^0 
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and, for the fluctuations. 
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where Vij denotes the second derivative with respect to 
Lpi and Lpj, and explicitly written as 
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We show two typical situations for the field evolu- 
tion. Figure |l| shows the result for the initial conditions, 
(/3i(0) = 10^ and (^2(0) = V^' where the prime denotes 
the derivative with respect to t. This corresponds to the 
circular orbit for the motion of the homogeneous mode. 
We can see that the instability band coincides exactly 
with that obtained analytically, since the upper bound 
of the instability band is 2/(^(0) = 0.002 in dimension- 
less units. On the other hand. Figure ^ shows the results 
for iy9i(0) = 10'^ and 1^92(0) = 0? which corresponds to the 
case that the homogeneous field is just oscillating along 
the radial direction. Almost all the modes are in the 
instability bands, but they have quasi-periodical struc- 
tures. This shows some features of the parametric reso- 
nance. Anyway, the most important consequence is the 
position of the most amplified mode, since it corresponds 
to the typical size of produced Q balls. Comparing with 
these two larger Q ball should form in the former 

(circular orbit) case, and the size will be about twice as 
large as that of the latter case, because the most ampli- 
fied mode kmax for the circular orbit case is about twice 
smaller than that for the just-oscillation case. 

We also see situations between circular orbit (e = 1) 
and just-oscillation (e = 0), where e represents the ellip- 
ticity of the orbit and defined by v?2(0) — \f2e. Higher in- 
stability bands become narrower and disappear as £ — *■ 1, 
and the lowest band also becomes narrower. 

Now let us move on to the fluctuations grown in the 
potential (||) where the gravity mediation effects domi- 
nate. Instabilities grow if K is negative, which is realized 
when the quantum corrections to the potential is domi- 
nated by the gaugino loops. We analyze this situation in 
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FIG. 1. Instability band for circular orbit e = 1 in the 
logarithmic potential (|l|). 
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FIG. 2. Instability bands for radial oscillation e = in 
the logarithmic potential (|l]). 

Ref. [|4j, so we only show the results here. The instabil- 
ity modes are in the range < {k/aY < 3\K\my^. The 

most amplified mode is (kmax/o.)'^ = 3|if|TO3^2/'^' ^^'^ 
the maximum growth rate is amax — 3|ii'|TO3/2/4 [|l^ . 

V. Q-BALL FORMATION 

Generally, Q balls are produced during the rotation 
of the AD field. In Ref. Q, we investigated the Q-ball 
formation in the gauge-mediation scenario on the three 
dimensional lattices, and found that it actually occurs. 
In there, we used the effective potential of the AD field 
as 

m=™^,log(^l + ^)+^|$r, (28) 

and took the initial conditions for the homogeneous mode 
as r(0) = 100, (pi{0) = A, ip[{0) = 0, ^2(0) = 0, and 



((32 (0) = B, where we varied A and B in some ranges, 
and added small random values representing fluctuations. 
However, at the initial amplitude which we took, the po- 
tential is dominated by the nonrenormalizable term, and 
for our choice of the initial time, the cosmic expansion 
is weaker than the realistic case. Thus, we could only 
find some restricted relations between initial values and 
charges of Q balls produced. 

Therefore, we use only the first term of Eg. (psf) , or 
equivalently, Eq.(|l^) with m{T) = const, for the poten- 
tial, and take initial conditions as 

(^2(0) = ^3, ip'M^S^ + h. 

(29) 

= 3^ = — 

where all variables are rescaled to be dimensionless pa- 
rameters, h = H /m, and e represents how circular the 
orbit of the field motion is e = 1 for the circular orbit 
and e = for the radial oscillation. (5's are fluctuations 
which originate from the quantum fluctuations of the AD 
field during inflation, and their amplitudes are estimated 
as 10~^ — 10"^" compared with corresponding homoge- 
neous modes. The initial time is taken to be at the time 
when H{t) — ui, which is the condition that the AD field 
starts its rotation. 

We calculate in one, two, and three dimensional lat- 
tices. In the first place, we set £ = 1, and vary the initial 
amplitude of the AD field ipo. It is usually considered 
that e may be very small in the gauge mediation sce- 
nario, because the gravitino mass 7713/2 is very small so 
that the amplitude of A-terms may be small. However, 
as we have shown above, e = 1 case can be achieved in 
certain situations. Figures ^ - |^ show the dependence 
of the maximum Q-ball charge produced upon the initial 
amplitude ipQ for one, two, and three dimensional lattices, 
respectively. 

We find the relation between the charge of the Q ball 
and the initial amplitude of the AD field: 

Qmax^ Pd^l'^'', d= 1,2,3, (30) 

where PdS are some numerical factors; Pi ~ 0.1, P2 ~ 
0.02, and /Sa « 6 x 10~"'. This relation can be understood 
as follows. Since the wavelength of the most amplified 
mode is as large as the horizon size, we can regard that 
one Q ball is created in the horizon size at a creation 
time. If we assume that Q balls arc created just after the 
AD field starts its rotation, the horizon size is H^^ ^ 
ijj^^ ^ (f)Q/ni^. Therefore, the charge of a Q ball is 

Q ^ H-^n^ ^ u-'^uj^l ^ m^-'' (31) 

This corresponds to the results of numerical calculations: 
Q oc fg'^'^- The prefactors Pd cannot be determined by 
the above analytical estimation. Actually, the formation 
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100 1000 10000 100000 1e+06 

FIG. 3. Q ~ ipo relation in the logarithmic po- 

tential on one-dimensional lattices. Dashed line denotes 

Qmax = 0.15<(9§. 



FIG. 5. Q — ipo relation in the logarithmic poten- 
tial on three-dimensional lattices. Dashed line denotes 

Qn^ax = 6 X W^ifiQ. 
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(p. 

FIG. 4. Q — ifio relation in the logarithmic po- 

tential on two-dimensional lattices. Dashed line denotes 
Qn,ax = 2 X IQ-Vo- 

time is a little bit later than the time when the rotation 
of the AD field starts, and the number of Q balls in the 
horizon size is more than one. Therefore, we will use 
/33 « 6 X lO""* for the estimation of Q-ball charge when 
we need later. 

Now let us see how the Q-ball formation occurs for 
e ^ 0. The charge of produced Q balls decreases with 
linear dependence on the charge density of the AD field, 
as shown in Ref. thus proportional to e. When the 
angular velocity of the AD field becomes small enough, 
both Q balls with positive and negative charges are pro- 
duced fl plJTsl . Therefore, the charge of the Q ball will 
be independent of small enough e. Using numerical calcu- 
lation, we find that the dependence of the Q-ball charge 
on e shows exactly the same features as expected. In 
Fig. ^, we plot the Q-ball charge on one dimensional 
lattices. We can see that Qmax is constant for small e 
where both positive and negative Q balls with the charges 
of the same order of magnitude with opposite signs are 
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FIG. 6. e-dependence of Q-ball charge in the logarithmic 
potential on one-dimensional lattices. The dotted line denotes 
Qmax — 1.3 X 10^£. On the other hand, the dashed line 
denotes the constant line Q ~ 3 x 10® . Closed squares denote 
the production of both positive and negative Q balls with 
charges of opposite signs and the same order of magnitude. 



produced, while linearly dependent on e around e ^ 1, 
where only positive Q balls are formed (smaller negative 
Q balls can be produced, but their charges are an order 
of magnitude or more smaller than the largest positive Q 
balls). We thus obtain the formula for the charge of the 
produced Q balls as 



Q 



ID 



3 X lo^ 

1.3 X lO^e, 



s < 0.2, 
e > 0.2, 



(32) 



in one dimension. 

Similarly, we plot the Q-ball charge on three dimen- 
sional lattices in Fig. 0. Also we can see that Qmax is 
constant for small e where both positive and negative Q 
balls are produced, while linearly dependent on e around 
e ~ 1, where only positive Q balls are formed. There- 
fore, we obtain the following formula for the charge of 
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FIG. 7. e-dependence of Q-ball charge in the logarithmic 
potential on three-dimensional lattices. The dotted line de- 
notes Qmax — 5 X 10*e. On the other hand, the dashed line 
denotes the constant line Q ~ 3 x lO'^. Closed squares denote 
the production of both positive and negative Q balls with 
charges of opposite signs and the same order of magnitude. 



the produced Q balls: 



Q 



3D 



3 X 10^ 

5 X lO^e, 



e < 0.06, 
e > 0.06, 



(33) 



The difference between the charges for e ~ 1 and e <C 1 
can be explained by the size of the produced Q balls. As 
seen in the previous section, the most amplified mode 
kmax of the fluctuations for the circular orbit {e — 1) is 
about twice as small as that of just-oscillation along the 
radial direction (e = 0). Therefore, twice larger Q balls 
are produced for e = 1 in one dimension, while about an 
order bigger ones are formed in three dimension, which 
exactly correspond to the difference between charges for 
£ ~ 1 and £ ^ 1. 

At very large amplitudes of the AD field, the gravity 
mediation effects for SUSY breaking will dominate, and 
the 'new' type of stable Q balls are produced As we 
see above, the potential is dominated by the form 



''3/2 



1 -l-i^log 



(34) 



In this case, the curvature of the potential does not de- 
pend on the amplitude so much, so the AD field starts its 
rotation when H ~ m.3/2. Thus, the initial time can be 
taken as t(0) = 2/3 if we rescale variables with respect 
to when we make them dimensionless. 

We have simulated the dynamics of the AD field on 
one, two, and three dimensional lattices, and find the for- 
mation of Q balls. Here, we take K = —0.01. (Actually, 
its value is estimated in the range from —0.01 to —0.1 
1^,^.) We plot the largest Q-ball charge in the func- 
tion of the initial amplitude (pn for one, two, and three 
dimensional lattices in Figs, y, and ^ respectively. 
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FIG. 8. Q — (fio relation in the potential (|34 

on one-dimensional lattices. Dashed line denotes 

Qmax = 0.2 X ipfj. 
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FIG. 9. Q — ifio relation in the potential (|34 

on two-dimensional lattices. Dashed line denotes 

QTnax = 0.05 X ipl. 



Similar to the logarithmic potential cases, we find the 
relation between the charge of the Q ball and the initial 
amplitude of the AD field as 



1,2,3, 



(35) 



where P^s are some numerical factors; /3i ~ 0.2, P2 ~ 
0.05, and (3^ w 0.006. This relation can be understood 
by the similar argument as the logarithmic potential case. 
The difference is that the time when the AD field starts 
oscillation and the most amplified mode (or instability 
band) enters the horizon do not coincide. Therefore, the 
charge of a Q ball is Il4| 



Q ^ H/m,:i/24>) - (|if|'/'m3/2)"''m3/2|i^|0^ 



"^3/2 



(36) 



where the subscript '/' denotes the values at the Q-ball 
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FIG. 10. Q ~ ifio relation in the potential (|3J) 

on three-dimensional lattices. Dashed line denotes 

Qn.a. = 0.006 X ifl 



formation time. This corresponds to the resuhs of nu- 
merical calculations: Q oc ipg. The prefactors I3d cannot 
be determined by the above analytical estimation. Actu- 
ally, the formation time is a little bit later than the time 
when the mode which will be mostly amplified enters the 
horizon, and the number of the Q balls in the horizon size 
is more than one. Therefore, we use /Ja « 6 x lO^'^ for 
the estimation of the Q-ball charge when we need later. 



VI. CHARGE EVAPORATION FROM Q BALLS 
IN HIGH TEMPERATURE 

A. Usual gauge-mediation type Q balls 

Since the Q-ball formation takes place nonadiabati- 
cally, (almost) all the charges are absorbed into produced 
Q balls l^jl^. Thus, the baryon number of the universe 
cannot be explained by the remaining charges outside Q 
balls in the form of the relic AD field. However, at finite 
temperature, the minimum of the free energy of the AD 
field is achieved for the situation that all the charges exist 
in the form of free particles in thermal plasma ||5|] . This 
situation is realized through charge evaporation from Q- 
ball surface 13|. In spite of the fact that the complete 
evaporation is the minimum of the free energy, the ac- 
tual universe is filled with the mixture of Q balls and 
surrounding free particles, since the evaporation rate be- 
comes smaller than the cosmic expansion rate at low tem- 
peratures. (As we mentioned in the Introduction, the free 
energy is minimized at the situation that some charges 
exist in the form of free particles (in thermal plasma) 
and the rest stays inside Q balls, if the Q-ball charge 
are large enough for the chemical equilibrium |p| . In this 
case, the charge of the Q ball should he Q ^ 7]g ^ 10"*".) 
The nonadiabatic creation of Q balls and the later charge 
evaporation takes place, since the time scale of evapora- 
tion is much longer than that of the Q-ball formation. 



Notice that the charge of the Q ball is conserved for the 
situation that the evaporation of charges is not effective. 
However, the energy of the Q ball decreases as the tem- 
perature of the universe decreases. This happens because 
the Q ball and surrounding plasma are in thermal equi- 
librium: i.e., the Q balls and thermal plasma have the 
same temperatures. Thus, the destruction process of the 
Q ball by the collision of the thermal particles, called the 
dissociation ||l^ cannot occur in actual situations. 

The rate of charge transfer from inside the Q ball to its 
outside is determined by the diffusion rate at high tem- 
perature and the evaporation rate at low temperature . 
When the difference between the chemical potentials of 
the plasma and the Q ball is small, chemical equilibrium 
is achieved and charges inside the Q ball cannot come out 
JlTf . Therefore, the charges in the 'atmosphere' of the Q 
ball should be taken away in order for further charge 
evaporation. This process is determined by the diffusion. 
The diffusion rate is given by ll^l 



dt 



-AnDRQfiQT^ An AT, 



(37) 



where D — A/T is the diffusion coefficient, and A = 4—6, 
fiQ ^ uj is the chemical potential of the Q ball. On the 
other hand, the evaporation rate is ||] 



dt 



f-^plasma 

J.2 



m{T) 



Q 



1/4 



(38) 



where fJ.piasma ^ fJ-Q is the chemical potential of thermal 
plasma, and iiq ~ cj ~ m{T)Q~^/'^ is used in the second 
line. m(T) and C change at T = as 



m(T) 



Therefore, we get 



evap 



dQ 
dt 



(T < m^), 

1 (r>m^). 



-47rTQi/^ (T > m^), 
2^4 



(39) 



-47r- 



(r<m^). 



(40) 



In order to see which rate is the bottle-neck for the 
process, let us take the ratio R = Tdiff/^evap- For T > 
TO^, the ratio becomes R = AQ~^/^. If i? < 1, the 
diffusion rate is the bottle-neck for the charge transfer, 
and this condition meets when the Q-ball charge is large 
enough as 



Q > 10^ 



(41) 



Since, as we will see shortly, we are interested in the 
dark matter Q ball, the charge should be as large as 
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10^^. More conservatively speaking, Q > 10^^ for stable 
against the decay into nucleons 1^. 

On the other hand, when T < m^, the condition R < 1 
corresponds to 



111 



(42) 



and the transition temperature is lower than for 
large enough Q-ball charge, which is again the interesting 
range for the dark matter Q balls. 

We must know the time dependence of the tempera- 
ture for estimating the total evaporated charge from the 
Q ball, which will be obtained by integrating Eqs.(|3^) 
and ( p^ ) [or (|4C|)]. Although thermal plasma (radiation) 
dominates the energy density of the universe only after 
reheating, it exists earlier, and this subdominant com- 
ponent will cause both the change of the shape of the 
effective potential for the AD field and the charge evap- 
oration from Q balls. The temperature of the radiation 
before reheating is 



1/4 



W/4 



(43) 



where F/ is the decay rate of the inflaton field, and the re- 
lation between the decay rate and reheating temperature, 
Trh ~ (-^r/)-^/^, is used in the last equality. Of course, 
when reheating occurs at ~ F/, the temperature be- 
comes as large as the reheating temperature: T ~ Trh- 
After reheating, the universe is dominated by radiation, 
and the temperature evolves as T oc H^/^ . 

Since t oc T~* before reheating, and t oc T^^ after 
reheating, the diffusion rate with respect to temperature 
T can be directly obtained from Eq.(p7[) as 



dQ 
dT 



10- ^" 



rp4 



diff 



10 



M 



(T > Trh).. 
(T < Trh). 



(44) 



On the other hand, the evaporation rate with respect to 
T is a little more complicated, and we can divide into four 
cases, depending on the temperature compared with the 
reheating temperature and the mass of the AD particle. 
The time dependence of T changes at T = Trh, and 
the rate dQ/dt changes at T — m^. Combining all the 
effects, we get 



lO^^Q^/^ {T>TRH,m,), 
lO^^Q'/' {TRH<T<m,), 



dQ 
dT 



10- 



T2 
MT 



(m^ < T < Trh), 



0"!^ iT<m^,TRH). 



(45) 



We can now estimate the total evaporated charge from 
the Q ball. According to the relation among T^,, Trh, 



and m^, there are three situations: (A) < < Trh, 
(B) n < Trh < (C) Trh < < m^. As we will 
see shortly, the evaporated charges in all cases are the 
same order of magnitude. Therefore, we will show only 
the case (A) here. For T > Trh, 



dT 



2^4 



(46) 



so we have 



AQ(T > Trh) = Q^nu - Qrh 



3MT 



RH \ rp3 rp3 

^ RH ^ init 



(47) 



where subscript 'init' denotes the initial values. If we 
assume Tinu ^ Trh, we obtain AQ{T > Trh) ^ 
3M/Trh. For T* < T < Trh, 



dQ M 

dT ^ r2' 



(48) 



so the evaporated charge in this period is estimated as 

(49) 



AQin <T< Trh) - lOM (7^-7^ 

\ J * J RH 

Finally, for T < T*, we must integrate equation 
^ ^ 10^QV4. 



dT 



We thus get 



AQ{T < T,) ^ 5^Qy^ 



(50) 



(51) 



where Qo is the amount of the Q-ball charge at present. 
Adding three evaporated charges, we have 

AQ-10— -7- — + 5 — —Qq. 52) 

T* Trh 

If Trh S> T^, the second term can be neglected. On the 
other hand, if Trh ^ T*, the sum of the first and second 
terms is the same order of magnitude as the first one. 
Therefore, from the first and second terms, we get 



-t* "10 



12 

init 



2.4 X 10' 



VTeV/ 



1024 



1/12 



(53) 



where we use Eq.(|4^) in the first line. For large enough 
charge such as Qmit S> 10^*, there is little difference 
between Qi„it and Qq, so the third term can be estimated 
as 

AQ(3) ^ ^^n^Q-^l^Q^ - S^OL^I^ (54) 



9 



which is the same order of magnitude as the sum of first 
and second terms. On the other hand, if Qmit < 10^^, all 
the charges are evaporated before the temperature drops 
to , and there is no contribution from the temperature 
below r*, which is the third term AQ^'^^ 

For the case (B), the reheating temperature is almost 
the same as the transition temperature T* , the amount of 



evaporated charge should be about the same as Eq.(53). 
Although it is not apparent in the case (C), we will now 
see that the amount of evaporated charge is as the same 
order of magnitude as in the case (A). Above the tem- 
perature T*, the diffusion rate determines the speed of 
the whole process, and it is exactly the same as the case 
(A). Thus AQ{T >n) ^ 3M/n. For Trh < T < T,, 
the rate is determined by the evaporation rate, and we 
should use 



dT ^ mlT 



)l/4 



(55) 



Integrating this equation, we obtain 



AQ{Trh <T<T,)^ io^Kq^ log (56) 



T, 



B. Gravity- mediation type Q balls 

Now we will show the evaporated charges for the 'new' 
type of stable Q ball The evaporation and diffusion 
rate have the same forms in terms of Q-ball parameters 
Rq and u!. The only differences are that we have to use 
the features for the 'gravity-mediation' type Q ball, such 



Q 



\K\ 



-1/2, 



'-3/2, 



"T-3/2, 



(61) 



diff 



be- 



and the transition temperature when T^vap 
comes = Ai/3|i4:|i/6(m3/2m2)i/3. As in the 'usual' 
type of Q balls, where the potential is dominated by the 
logarithmic term, the charge evaporation near is dom- 
inant, and the total evaporated charges are found to be 



AQ 



10- f ^) 



-1/3 



VTcV/ 



-2/3 



(62) 



RH 



VII. COSMOLOGICAL Q-BALL SCENARIO 



On the other hand, for T < Trh, we must integrate 



dT 



and the result is 



AQ{T < Trh) - 10 



MT 



RHqI/A 



(57) 



(58) 



Since Trh ^ T*, Equations ( |56| ) and (|5^) have almost 
the same form, we can regard the charge evaporation in 
these stages as 



AQ - 10 



MT 



RHqI/4 



(59) 



For the contribution from this term to be dominant, we 
must have condition Trh ^ mrpQ~^/^^ . However, since 
Trh ^ T,: ^ TO^Q~^/^^, contributions from this stage 
is at most comparable to that in the region T > T^. 
Therefore, we can adopt the estimation of the evaporated 
charge from the Q ball as 



VTeV/ VlO 



AQ ^ lO^Qi/- ^ 2.4 X 10i« ( A 



for any cases. 



1/12 



(60) 



A. Usual gauge-mediation type Q balls 

Now we would like to see whether there is any consis- 
tent cosmological scenario for the baryogenesis and the 
dark matter of the universe, provided by large Q balls. 
In the first place, we will look for the situation in which 
the logarithmic term dominates the AD potential, and 
the 'usual' gauge-mediation type of Q balls are formed. 

Speaking very loosely, we know that the amount of 
the baryons in the universe is as large as that of the dark 
matter (within a few orders of magnitude) . In the Q-ball 
scenario, the baryon number of the universe should be 
explained by the amount of the charge evaporated from 
Q balls, AQ, and the survived Q balls become the dark 
matter. If we assume that Q balls do not exceed the 
critical density of the universe, i.e., fJg ^ 1, and the 
baryon-to-photon ratio as ?/b ~ 10^^*^, the condition can 
be written as 



riB 



enqAQ epqAQ epcfiflgAQ 



UjMq 



(63) 



where ilg is the density parameter for the Q ball. Using 

3.3x 



)3/4 



Pc.o 8h^ X 10-*^GeV*, and n 



7,0 



Mq ~ . 

10^'^^GeV'^, where /io(~ 0.7) is the Hubble parameter 
normalized with lOOkm/sec/Mpc, we obtain the ratio of 
the charges evaporated to be the baryons in the universe 
and remaining in the dark matter Q ball: 



AQ 



10"e-i77B^ 



1/4 
-1 / "10 \ 

? VTevJ 



)-l/4 



(64) 
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As we mentioned earlier, in the inflaton-oscillation 
dominated universe between inflation and reheating, the 
temperature has different dependence on the time (or the 
Hubble parameter) from that of the radiation dominated 
universe, and it reads as, from Eq.(B3), 



{MTf 



RH 



(65) 



At the beginning of the AD field rotatioUj H ^ 

m'^{T)/(j)o ~ rV0o- Inserting this into Eq.(|5|) and 
rephrasing it, we obtain 



(66) 



In the Affleck-Dine mechanism with Q-ball production, 
the baryon number of the universe can be estimated as 




VB 



nB,RH 



Pimh/Trh Pij/Trh 



(67) 



where subscript RH and / denote the values at the time 
of the reheating and the formation of Q balls, respec- 
tively. Notice that ub and pi are proportional to a^^ 
during the inflaton-oscillation dominated universe before 
reheating. Thus we must have the baryon number den- 
sity at the formation time which will be evaporated later. 
It can be written as 

riBj' ^ TBen^j ~ VBSujcj)^ 

^ lO^'^Bn-Q^ (^) TVoQ-^/^ (68) 

where we use Eq.(|^, and uj ^ m^(T)/0o ^ T'^ /4>o- On 
the other hand, the energy density of the inflaton when 
the AD field starts its rotation is 



9o 



(69) 



Thus, from Eq.(|67|), we can write the baryon-to-photon 
ratio as 



VB 



Q VTcV/ TrhM^ 



-1/4 



(70) 



where Eqs.(|66|), (p8[), and ( |69| ) are used. We thus get the 
Q-ball charge as the function of the initial amplitude of 
the AD field and the reheating temperature as 



44 fm^V <t>l^ 0-4 



(71) 



Now we will obtain another relation between the initial 
amplitude of the AD field and the reheating temperature. 
From the analytical and numerical estimation for the Q- 
ball charge, we have 



T 



(72) 



for £ ~ 1. For e ^ 1, we have to replace (3 hy f3' = 7/3 
with 7 ^ 0.1. 

Since we obtain two different expressions for the Q- 
ball charge, we can get the initial amplitude of the AD 
field from Eqs.(^l]) and (72) as, taking into account the 
e-dependence when e ~ 1 in Eq.(F3), 



4.6 X lO'^e 



13^1/10^2/5 



Q 



_JL_)"'°(JIl±y"G.V. (73) 

6 X 10-4 / VTcV/ ^ ' 



Inserting this equation into Eq.(|6q), we get the corre- 
sponding temperature when the AD field starts the rota- 
tion. It reads as 

T ^ 2.3 X lO^e-i/^Of]-^/^ ( -^^) 

/ o \ -1/20 , /r 

xf ^ f^)''^GeV. (74) 

V6xl0-V VTeV/ ^ ' 

We also obtain the charge of the Q ball 
Q 9.3 xiU e l^^oSGevj 



6 X 10- 



Tim-""- <™) 



where we use the numerical estimation ([74), and insert 
Eqs.(|66|) and ( [73|) into it. Notice that the charge Q does 
not depend on the amount of the baryons, so that it just 
represents the charge of the dark matter Q ball. 

We must know how circular the orbit of the AD field 
motion is. It is necessary to obtain the expression for e in 
terms of other variables, say, the reheating temperature 
TnH ■ In addition to the condition of the amount of the 
evaporated charge for explaining both the baryons and 
the dark matter simultaneously, which can be seen in 
Eq. (|6^ , we also have the estimation of the evaporated 
charge from the Q ball (pOh. Equating these two, we have 



1.2xl0V/V/^<f^) 



-3 



Q \TeVJ 
Compared with Eq.(ra), e should be 



(76) 



1.3 X io"r2 



Q 



VB 

10-1' 



-15 



VlO^GeVy 



f3 



6 X 10-4 



VTeV/ 



(77) 



for explaining the amount of the baryons and the dark 
matter of the universe simultaneously. In the case of 
e < 1, Eq.(^ has no e-dependence, so Eqs.(^) - (|7|) 
also have no e-dependence. Thus, we instead obtain the 
formula for e as 
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£ ~ 1.5 X 10-^fl^r{^ 



/3 



Vio^Gevy 

16/15 



6 X 10-5 



8/3 



TeV 



2/5 



(78) 



To summarize, the conditions of parameters for the Q- 
ball baryogenesis and dark matter scenario to work are 
given as 

<^o ~ 4.6 X 10iV/i°f^|/' 

.(^)'""(^)-CeV, ,™, 



^6 X 10-4 
T- 2.3 X 10^e"^/2"l7-^/^ 



TeV. 

Trh \ 
lO^GeV / 



/3 



6 X 10- 



-1/20 



VTeV/ 



y.dxlU e Uq lioSGevi 



6 X 10- 



8/5 



VTeV/ 



-12/5 



1.3 X lO^^fln^ , 

^ VlO 



-Y 

Vio^Gevy 



6 X 10- 



-16 



/ TO^y 
VTeV/ 



(80) 



(81) 



(82) 



where we should omit e-dependences in Eqs.([79|) — <^^) 
for £ < 1. 



As we mentioned in Sec. [11, the initial amplitude of the 
AD field is determined by the balance between the Hub- 
ble mass term and the nonrenormalizable term. When 
the AD field starts rolling down its potential, the am- 
plitude becomes 0o ~ {HoscM^~^y/^"~'^^ for the super- 
potential W ~ 0"/M"-3, where Hose ~ Using 
Eq.(n6|), we can write it as 



00 - (r2^M"-2)i/« 



(83) 



We will see the range of n in which we can obtain the 
consistent scenario naturally. At the first place, let us 
consider the case n = 4, where e <C 1 will be led below. 
In this case, the initial amplitude becomes 



(84) 



Therefore, the required reheating temperature should be 



Trh 8.8 X lO^r^i''^ 



(3 



6 X 10- 



1/5 



VTeV/ 



-4/5 



GeV, 
(85) 



T - 2.0 X 10"f7!V^ 



6 X 10- 



\ 3/20 ,/r 



(86) 



e~8.8xl0-«17o''^^^f^ 



p 



6 X 10- 



8/15 . ^ , . 38/15 



VTeV/ 



so that all the charges evaporate and no dark matter Q 
balls exists. For the n — 5 case, 

and the reheating temperature should be 



6x10-4/ VTcV/ 



(90) 



which leads to, taking into account the fact that e ^ 1 
in this case, 



T - 8.3 X lO'nj/' ( ^ 



^ V 6 X 10-4 



Q - 5.5 X 10^^f7 



-8/5 / P 



3/5 



8/5 



(92) 



^6xl0-4y VTeV/ 
s^l.3xlO-V^/-(l^) 

(93) 



6 X 10- 



TeV 



In order for the Q ball to survive from the evaporation, 
the initial charge of the Q ball should be large enough. 
This condition is Q ^ AQ, and can be achieved from 
Eq.(|0l) if 



Q>7.4xl0"f^) 
VTcV/ 



-12/11 



(94) 



Imposing this condition on the required Q-ball charge 
above, we obtain the required mass of the AD particle as 



> 6.2 X 10-^17 



3o22/37 / P 



^ V6xl0-4 



-33/148 



GeV, (95) 



so that the degree of the ellipticity of the orbit of the AD 
field motion should be 



e > 3.5 X lO^^n, 



' Vio-i"/ 



P 



6 X 10- 



-21/74 



(96) 



where Eq. ([79[) is used. Then, for this reheating temper- 
ature, we have 



Since, e ^ 1 in this case, we have to use an order of 
magnitude smaller value for (3. Then, we get the following 
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values for the parameters in order for the Q-ball scenario 
to work naturally: 

~ 1.0 X 104/3,-33/148^^22/37 g^^^ 

e~6.7xl0-V40/3;''/''V''^"\ 
T ~ 8.3 X lO^/Sl^/^^^f^Q GeV, 
Tb,h ~ 2.0 X 10"^ [3f'^^n^^''^'^ GeV, 
,^0 ~ 1.5 X 10i3/3j/37j7^/37 Gev^ 
Q ^ 5.5 X 1016/39/37^^^24/37^ ^97) 

where (3^ = (3/ {6 x 10"^) and 7jb,io = vb/10^^°- Ac- 
tually, this parameter set is the lower limit for e, and 
larger values are also allowed. As we see later, the up- 
per bound comes from the condition that the Q ball 
is stable against the decay into nucleons. The allowed 
range is 6.7 x 10-3 < e < 5.1 x 10-^, or, equivalently, 
1.0 X lO-^GeV < < 2.0 x lO" GeV. 

Now let us move on to the n = 6 case. Repeating 
the similar argument, we obtain the consistent values for 
parameters 



e ~ 0.9717q39/11^15/11«4/11_54/11 



1 15/11^4/11 54/11 
''^B,10 Pi "^-0,430' 



T ^ 9.1 X 106ei/4f7Q^i/4^-i^^ Qgy^ 

TnH - 4.7 X 10^3/10^^^/5^3/10^-6/5 g^y^ 
00 - 6.4 X 10l3£VlOf^^/5^1/10^-2/5 g^y^ 

g ^ 1.5 X l^^'e^''nr''PT<L. 



(98) 



where Pi = /3/(6x lO"'^) and m0,43O = "1^/(430 GeV). In 
this case, the AD field rotates circularly and produce the 
baryon number maximally. In general, however, e <C 1 
cases are also allowed. Numerical calculations reveal that 
the m^-dependence of the Q-ball charge changes around 
e ~ 0.1 using the following results shown in the previous 
section: Q is proportional to e for e > 0.1, while constant 
for smaller values, which reflects the production of both 
positive and negative charged Q balls. Therefore, the 
charge can be written as 



Q oc 



C'" (e>0.1). 



(99) 



and the lower bound may be determined by the possible 

lowest mass, ^ 100 GeV. 

For the n = 7 case, we need extremely huge e such as 
6 X 10^, which cannot be realized. Higher n makes the 

situation worse. 

In order for the scenario to work naturally, we must 

check that the values of e are consistent with the A-terms 



derived from the nonrenormalizable superpotential. As 
we derived in Sec. Ill, the formula for the e is written as 



^3/200 
2^2 



(100) 



For n = 5 case, putting the values for 0o and T in 
Eg. (p7|) , we have 



e - 2.2 X 10-3f7 



/ o \ -21/37 
-30-18/37 ( P \ ( ^3/2 

6 X 10-5 / V GeV 



(101) 



Comparing this with the value of e in Eq.(|97|), we get 



^3/2 



0.33a 



o \ -21/74 ^ 
-479/111 I P \ f VB 



6 X 10- 



10 



-10 



GeV. 



(102) 



Therefore, we get a consistent Q-ball scenario naturally 
in the 7i = 5 case. Notice that the reheating temperature 
is low enough for m3/2 ^ 1 GeV to avoid the cosmolog- 
ical gravitino problem Notice that, for the larger s 
cases, 1713/2 gets too large in the framework of the gauge- 
mediated SUSY breaking, and the allowed range becomes 
~ 10 - 14 TeV. 
In the n = 6 case, we have 



o.83r2; 



-8/7 



6 X 10-4 



-2/7 



/ m0 \8/7 / TO3/2 \5/7 

U30GeV/ VGeV/ 



(103) 



for the values of 0o and T in Eq.(|98D, and comparing 
with the value of e in Eq. (|9^) , we get 



nT-3/2 



1.2 n: 



-51/21 



6 X 10- 



10/11 



/ X 297/55 / \ 21/11 

UsOGeV/ V 10-10/ 



GeV. (104) 



Therefore, we obtain the consistent scenario for n — 6, 
which is achieved, for example, if we choose the udd flat 
direction for the AD field and use W ~ (udd)^. Notice 
again that no gravitino problem exist in this scenario. 
The scenario also works for smaller e for 7723/2 ~ 1 GeV, 
and the allowed range of the parameter is 2.8 x 10-3 < 
e < 1, or, equivalently, lOOGeV < 771^ < 430 GeV. 

Now we must see whether the effective potential of the 
AD field is really dominated by the logarithmic term over 
the gravity-mediation term at this field amplitude. This 
condition holds if 



> 



'-3/2' 



io. 



(105) 



or, equivalently, 
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T > (m3/2</>o)'/' 



log 



y2 



-1/4 



(106) 



If we use the values in Eqs.(|9^) and (98), the right hand 
sides becomes ~ 1.0 x 10^ GeV and - 3.8 x 10^ GeV, 
respectively. Thus, the condition is met, and we have 
consistent cosmological scenarios in the gauge-mediated 
SUSY breaking model for the effective potential domi- 
nated by the (thermal) logarithmic term. 



B. Gravity- mediation type Q balls 

For the thorough investigation, we should consider 
whether the 'new' type Q-ball scenario works for low 
enough reheating temperature avoiding the thermal ef- 
fects. We will follow the same argument as we did for the 
'usual' gauge-mediation type Q ball. From the baryon- 
to-photon ratio (|6^), the ratio of the charge evaporated 
from the Q ball to that remained in the Q ball can be 
estimated as 



AQ 



where we have used Mq ~ TO3/2 Q an d put e = 1, since it 
is a natural realization (see Sec. III). Since the baryon- 
to-photon ratio can also be written as 



nB,osc 



Pl.osc 



RH 



m 



A'P/Ti 



3/2 



RH 



the reheating temperature can be estimated as 



Trh ~ 10-^ n 



Q 



M 



GeV. 



(108) 



(109) 



In the previous section, we have obtained the relation 
between the formed Q-ball charge and the initial ampli- 
tude of the AD field as 



"13/2 



(110) 



where /3 « 6 x 10^"^ from our simulations. On the other 
hand, we have the condition that the evaporated charge 
and survived stable Q balls explain for the baryon and 
the dark matter of the universe, respectively, as |^ 



VMeVy 



VTeV/ 



(111) 



Therefore, from Eqs.( |llQ| ) and (111), we obtain the re- 
quired amplitude of the AD field as 



do ~ 4 X lo'^ny' ( ^ 




-1/2 



^m3/2_Y''^ / N -1/3 



Inserting this into Eq.( |l09| ), the required reheating tem- 
perature is 



Trh - 3.5 X 10 



7 P 



6 X 10-3 



When the gravity-mediation term rriy^'t''^ dominates 
the AD potential, the field starts its rotation aX H ^ 
™3/2j so the temperature at that time is Tqsc ~ 
(AfTl^mg/a)^/''. Thus we have 



Tosc - 4 X 10 



7 P 



1/2 



6 X 10-3 



^(Si)"""(^)"<=- <'»' 

We must obtain the constraint that the gravity-mediation 
term dominates over the thermal logarithmic term in or- 
der for the scenario of the 'new' type Q ball to be suc- 
cessful. It reads as 

Ti,\og(^) <ml/,cl,l (115) 

\ osc / 



and we can rephrase it as 



1/3 



(116) 

Therefore, this range of the gravitino mass is too large for 
the gauge-mediated SUSY breaking model, and the ther- 
mal logarithmic term dominates the potential at these 
field values. We thus find no consistent model for the 
'new' type of stable Q balls. 



VIII. GENERIC MODELS FOR GAUGE 
MEDIATION 

A. Dominated by zero-temperature potential 

As we mentioned earlier, the scale of the logarithmic 
potential could be much larger than in the general 
context of the gauge-mediated SUSY breaking. For the 
complete consideration, we will discuss this situation in 
this section. Here we will express it as 
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V 



Mi 



'> Ms), 
'< Ms), 



On the other hand, when T < T^,, the charge transfer is 
(117) determmed by the evaporation rate, and the evaporated 
charge can be estimated as 



where Ms is the messenger mass scale. Formation of the 
dark matter Q balls will take place at large field ampli- 
tudes, so the mass, the size, etc. have to be different, 
and their expressions are as follows: 



M, 



Q 



OJ ' 



M| 



(118) 



In order for this potential to dominate over the thermal 
logarithmic potential, we need condition Mp > T. Oth- 
erwise, the results discussed in the next subsections have 
to be appHed. Since this type of the Q ball should be sta- 
ble against the decay into nucleons, the Q-ball mass per 
unit charge must be smaller than 1 GeV. This condition 
holds for 



Q>io24 ( Mp y 



(119) 



We can regard that all variables are rescaled with re- 
spect to A'lp in our simulations now. So the charge of 
the produced Q ball is 



Q = 



Mf 



(120) 



where /3 « 6 x 10 * again. The evaporation rate has to 
be changed to 



Uvap- 



-47r— (T > m^), 

% (121) 



Depending on the Q-ball charge, the transition tempera- 
ture T* where the diffusion and the evaporation rates are 
equal can be written as 

r amfQ-^/^ in > m^), 

* \ Ai/3(^2^^)i/3Q-i/i2 (T, <m^), 

where A ~ 4 is defined in Eg. (p7[) . These two temper- 
atures coincide when Q = Qcr ~ A^Mpm'^'^ . Taking 
the stability condition ( |119| ) into account, we must con- 
sider only the case with Q > Qcr, which corresponds to 
r» < 777.0 case. 

The charge transfer from inside the Q ball to its outside 
is determined by the diffusion rate when T > T*, so the 
evaporated charge during this period is 



Ag(T > T,) ^ IQ^ 

1 * 



4.6 X 10 



15 / ™0 Y"^'^ ( 



Mp 



-1/3 



AQ(T < T,) ^ -^Q^'\ 



4.6 X 10' 



102GeV 



VlO^Gev) 



-1/3 



)1/12 



(124) 



Therefore, the total charge evaporated from the Q ball 
can be written as 



AQ ~ 4.6 X 10 



15^ f Mf V'^' 1/12 

102GeV/ ^lOeGeVy ^ ' 



(125) 



Now we can impose the survival condition, which im- 
plies that the Q ball should survive from the charge evap- 
oration and become the dark matter of the universe. It 
reads as Q ^ AQ, and equivalent to 



> 1.2 X 10'^ 



\-8/ii / Mf 
102GcV/ IvlOSGeVy 



7770 



(126) 



Since the baryon number and the amount of the dark 
matter is related as in Eq.(|6^), we get the charge of the 
Q balls for this type as 



Q^3.2xl0X'V/^(l|^'"'' 



/ m^, \-i / Mf \ 
UO^GeV/ VlO^GcV/ 



(127) 



We call this the baryon-dark matter (B-DM) condition. 

In addition to three constrain ts m entioned above, there 
is another one. The potential (117) at large field values 
must dominate over the thermal logarithmic potential, 
and this condition is set by Mp ^ T, as we mentioned 
above. Rewriting it, we obtain 



Mf > 2.7 X lO^e 



7^-5/12^^1/36 



f3 



7/18 



6 X 10-4 

5/12 / 777a \ 5/18 



where we have used the required initial amplitude of the 
AD field. 



00 ~ 1.5 X 10"e3/8f^i/36 ( P \ 



-1/4 



10-1' 



-3/8 



I lOSGeV J ■ ^^^^^ 

which is derived by equating Eqs.( [l20|) and (|l27|) 
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FIG. 11. Constraints on {Q,Mf) plain with = 100 
GeV. There are four conditions: (a) non-dominant thermal 
logarithmic potential, (b) stability against the decay into nu- 
cleons, (c) survival from charge evaporation, and (d) right 
relation between the amounts of the baryon and the dark 
matter. 



We can now put these four constraints on the param- 
eter space {Q, Mp) as shown in Fig. |ll[ In this figure, 
there are four constraints: line (a ) is the condition that 
the potential is dominated by Eq.(117), and the thermal 
logarithmic potential is negligible, line (b) is the stability 
condition where the Q ball cannot decay into nucleons 
and becomes the dark matter of the universe, line (c) 
represents the survival condition that the Q ball should 
survive from charge evaporation, and line (d) denotes 
the B-DM condition that have the right relation between 
the amounts of the baryons and the dark matter. We 
also show arrows to notice that which side of these line 
are allowed. As can be seen, there is no allowed region in 
the parameter space if we combine four conditions above, 
which implies that this type of the Q ball is very difhcult 
to explain both the amount of the baryons and the dark 
matter simultaneously. 



B. Dominated by finite-temperature potential 

Here we will investigate for the Q-ball formation in 
the generic potential, but the thermal logarithmic term 
dominates the potential at the formation time. We have 
considered the special case Mp = in the previous sec- 



tion, and the results which we have derived above can be 
used if we replace by Mp appropriately. As we men- 
tioned in the last subsection, the only exception where 
the simple replacing cannot be ap plied is the estimation 
of the evaporated charges [see Eq.( |125| )]. Similar analysis 
reveals that the formation of the Q balls takes place with 
the following parameters: 



00 - 2.9 X 1012^1/10^^2/5^1/10^^-2/5 g^^^ 



(130) 



T - 9.2 X 10'£"i/2°f7^i/^/37^^^°A45TflH^5 GeV, (131) 



Q ^ 5.9 X lO'^n'^^'p^^'M-''/' 



rp — i 

F,6 ^i?,ff,5' 



(132) 



e ^ 2.2 X 10''n-^Pi''vB'>^'2Mk6TkV5^ (133) 

where Pi = /3/(6 x IQ-"'), Mpe = Mf/{10^ GeV), 
Trh,5 = Trh/{10^ GeV), t]b,io = r/s/lQ-i", and m^,2 = 
TO^/(102GeV). 

The initial amplitude of the field is determined by the 
balance between the Hubble mass and nonrcnormalizable 
terms. Since the scenario naturally works only for rt = 6 
case, we depict the results for this case only. Since the 
initial value of the field is 

^0 - [Tl^M^f"" ~ 8.3 X lO^^T^/^^^ GeV, (134) 

we get the relation between the reheating temperature 
and Mp as 



r«H,5 - 4.3 X lQ-^e^l^'n%'(3l"%I-'^\ 
where we have used Eq.( |l30| ). Then we have 

T - 4.0 X IQ^e^'^nqPl'^Mpl GeV, 

Q ^ 1.7 X 103V/5^]-l'/%2/^M^^/^ 
s ^ 2.5 X 10-f7-^/"/3Vii^i5/n^,^^^^^^ 



(135) 

(136) 
(137) 
(138) 



For e <C 1 we have Eqs.(136) and (137) without e- 
dependence, and 

e ^ 3.6 X 10917o"/^/34/i5,,^_^„mWi5 2/3^ ^^39) 



W e get the right answer if we equate Eqs.(127) and 
(137), and the required values are as follows: 



Q ~ 1.6 X 10^ 



Mf 



4.5 X 10-^GeV, 



1.0. (140) 



As is the case of the specific potential discussed in the 
previous section, the Q-ball charge depends linearly on 
e for e > 0.1, and constant for e < 0.1. Therefore, we 
obtain the iW^F-dependence of Q-ball charge as 



M 



12/5 



Mp > 2.5 X 102 GeV, 
Mp < 2.5 X 102 GeV, 



(141) 



in the range Mp = 10^-4.5 X 102 GeV, or, equivalently, 
£ = 6.7 X 10^3 _ i Q foj. = 100 GeV. Notice that the 
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condition T > Mp must be hold for this scenario to work 
properly. It can be rephrased as 



Mf < 6.2 X GeV, 



(142) 



where Eqs.(131) and (135) are used. Therefore, the above 
scenario for n = 6 is allowed by this condition. 



IX. DELAYED Q-BALL FORMATION 

Since we are interested in the Q-ball formation, we 
have considered only K < Q cases in the gravity-mediated 
SUSY breaking effects. However, it is possible for some 
flat directions that K become positive. As opposed to 
the gravity-mediation model, the sign of the K term has 
an ambiguity because of the complexity of the particle 
contents in the gauge-mediated SUSY breaking model. 
Here, we consider this possibility and will find some con- 
sistent scenario for the dark matter Q balls providing the 
baryons of the universe. We will derive constraints in the 
generic potential first, and put Mp = for the specific 
one later. 



A. Dominated by the thermal logarithmic potential 

Let us investigate the case which the gravity-mediation 
and thermal logarithmic terms dominate the effective po- 
tential at large and small scales, respectively. The critical 
amplitude of the AD field is determined by 



''3/2V 



Thus, we get 



T2 

eg 



(143) 



(144) 



where the subscript 'eq' denotes the variables which 
are estimated at the equality of two different potentials 
above. At this time the horizon size becomes larger than 
that at the beginning of the field rotation. At larger 
scales when the gravity-mediation term dominates, the 
amplitude of the field decreases as oc a~^^^ as the uni- 
verse expands. The universe is assumed to be dominated 
by the infiaton-oscillation energy, a cx t^^^ oc H~^/^, 
which leads to 6 (x H. We can thus find the horizon size 



aX 6 — 



^eq ■ 



It reads as 




(146) 



Since the typical size of the resonance mode in the in- 



stability band is k, 



1 

res 



eq 



''3/2- 



the horizon 



size is larger, so that the instabilities develop as soon 
as the field feels the negative pressure due to the log- 
arithmic potential, and Q balls are produced with the 
typical size ~ k~^^. Since the angular velocity is written 
as LUeq ^ {T'^ql 4>eq)~^ ~ ^3/2 at this time, the number 
density of the field can be expressed as 



^eq4>lq 



rp4 

eq 
"^3/2 ' 



(147) 



This corresponds exactly to the number density just di- 
luted by the cosmic expansion: 



He 



an 



rp4 

eq 



(148) 



'■eq / \ V^eg / W3/2 

Therefore, the charge of the produced Q ball will be 



eg 



"^3/2 



"^3/2 



eg 



(149) 



where the last and the second last terms have the same 
forms as the formulas of the charge estimation for the 
'usual' and 'new' type of the Q balls, respectively. We 
thus apply the numerical results which we obtained for 
the 'usual' type, and assume the charge of the Q ball as 



T, 



eq 



(150) 



where /3 « 6 x 10^^. ^ Notice that e ~ 1 in this case, 
since the oscillation (rotation) of the field starts when 
the potential is dominated by the gravity-mediation term 
where u) W3/2, as mentioned at the end of Sect. III. 



The temperature of the radiation at 



mined by 



eq 



^eq 



is deter- 



We thus find it as 




RH- 



(151) 



The baryon-to-photon ratio can again be expressed in 
two ways. One is 



eq 



Ho 



^eq 



(145) 



where subscript 'osc' denotes the values at the beginning 
of the field oscillation (rotation) and Hose — 'mz/2- We 
find the horizon size &t 6 = Sea as 



may be larger than this value, since the cosmic expansion 
is weaker than the situations which are done in numerical cal- 
culations above, and the Q-ball formation takes place earlier. 
However, the following estimates do not change because of 
the weak dependence on 
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and we can rewrite the ratio rs as 
AO ^„14 o-i f 

On the other hand, we have 



Q 



-1/4 



Pl,osc 



rB4>lTRH 



RH 



m3/2 



(152) 



(153) 



(154) 



Inserting Eqs.( |l50D , ( [l5l|) , and ( 153 ) into this equation, 
we obtain the initial amplitude of the field: 



2.9 X 10i2f]|/5/3yio 



M 



-2/5 
F,6 



GeV. 



(155) 



Then the temperature of the radiation and the amplitude 
of the field at cj) — (t)eq are 



Te, ^ 9.2 X 10^f^o'/'/37'/'°A4/'r;^ff,5 GeV 



,.4 X iQi^riQ^/^/?,-^/^" 

2/5. 



F,6^RH,5 



"^3/2,GcV GeV, 



(156) 



(157) 



respectively, where ?ti3/2,goV = "^3/2/GeV. We can thus 
estimate the charge of the Q ball as 

g ~ 4.3 X lO'^n'^'p^'m^'^T'^H^.m-;^^^^^ (158) 



which is 



(159) 



Notice that we need the constraint 00 ^ 
rephrased as 

00 > 1.5 X 10'^TRH,rom-^^^^^^ GeV 

for this situation to take place. 

There are two other conditions to be imposed. One of 
them is that the temperature of the radiation at = c/jeq 
is larger than Mp. Otherwise, the thermal logarithmic 
term in the potential is not the dominant one, which 
contradicts the assumption of this subsection. This con- 
strains the reheating temperature from below, such as 

Trh > 1.1 X 10^f7o''^^y^°M^/gGeV. 



(160) 

The second one is that the gravitino mass should be 
smaller than 1 GeV, since we are discussing the gauge- 
mediated SUSY breaking model. As will be seen shortly, 
this constrains the allowed region of Mp from below, 
when we assume that the initial amplitude of the AD 
field is determined by the balance between the (nega- 
tive) Hubble mass and nonrenormalizable terms. From 
this assumption, we get 0o ~ (m3/2M"^'^)-'^/("^^^ For 
some n, we can write them down as 



1.5 X lO^rriy^ GeV {n 

1.8 X 10i2^3^2.GcV GeV {n 

6.1 X 10^7713/2 GeV {n 

5.1 X 10"m3^2 GeV {n 

2.1 X 1015to3^2,GoV GeV (n 



4) , 

5) , 

6) , 

7) , 



(161) 



Equating these with Eq.(155) and using the condition 
that the gravitino mass should be less than 1 GeV, we 
obtain the constraints on Mp. 

Combining these with the condition that the reheating 
temperature must be larger than ^ 10 MeV for the sake 
of the successful nucleosynthesis ||20| , we obtain the lower 
limit of Mp for the allowed range, using Eq.(160). In 
general, parameter space (Q, Mp) is constrained by three 
conditions: the stability, survival, and B-DM con di tions , 
which we showed in the last section. See Eqs.( |ll9] ), ( 126 ), 
and (127). We find the allowed range of Mp only for 
n = 6 and n = 7, if we take into account all the conditions 
which we have mentioned above. These are 

Mi. - 4.9 X 10^ - 1.1 X 10* GeV (n = 6), (162) 
Mf - 1.0 X 10^ - 2.5 X 10^ GeV (n = 7). (163) 

For n — 6, the lower and upper limits come from the 
conditions 7713/2 ^ 1 GeV and Mp < Teq, respectively. 
On the other hand, for n = 7, the upper limit comes 
from the condition Trh 10 MeV, and we assume that 
Mp > 100 GeV. Therefore, we have consistent scenarios 
in this model, and, for the above ranges for Mp, the 
allowed parameter regions are 



Q - 1.3 X lO^'' - 2.8 X 10^\ 
Trh - 1.0 X 10^ - 29 GeV, 



for 71 = 6, and 
Trh ' 



3.2 X 10^^ - 5.1 X 10^^ 
62 GeV - 10 MeV, 



(164) 
(165) 



(166) 
(167) 



for n = 7, if we use the B-DM condition in order for 
the dark matter Q balls to supply the baryons in the 
universe. The allowed regions will be plotted in Figs. |l^ 
and |l^ in the next section, where the constraints from 
several experiments are also plotted. 



B. Dominated by zero temperature generic 
logarithmic term 

Now we will consider the situation when the tempera- 
ture of the radiation is rather low, and the effective po- 
tential is dominated by the generic logarithmic term 



Y 



M| log 



Ml 



(168) 



Similar discussion follows along the line which we made 
in the last subsection. The critical amplitude of the field 
is determined by 




ZjlVeqi 



(169) 



and we thus get 0eg ^ M'j^/m3/2- 
is then written as 



The horizon size at 



^eq 
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H ^ Ml 



(170) 



charge of the produced Q ball [Eg. (171)] should satisfy 
the stabihty condition [Eq. (|ll9| )] , as 



which is larger than the typical resonance scale in the 
instability band: k^J^ ^ {Mp/(j)eq)~^- Therefore, the 
field feels spatial instabilities just after its amplitude gets 
smaller than (/)eq. The charge of the Q ball can be esti- 
mated as 



^ ^ ^ lO^O/J.M^ ^mg/, Q^^. (171) 



^3/2 < 0.16/5^^ GeV. (178) 
Therefore, we obtain the allowed values of parameters as 



m3/2 < 0.16 GeV, 

10 MeV < Trh < 110 GeV, 

Mf > 100 GeV, 



(179) 
(180) 
(181) 



Using Eqs.(|l53|), (|154|), and ( |l7l| ), we obtain the ini- 
tial amplitude required for the Q-ball formation in this 
scenario as 



^0 9.5 X 10'"r!j/'/3y'r^^{6 GeV, 



(172) 



where Trh,6 = Trh /{W'^GeY). 

At the time of the production of Q balls when (jj ~ (peq, 
the temperature of the radiation can be estimated from 

Teq - {MTRHHeqy/'^, SO that we have 



Te, - 7.1 X 10'r!Q'/'/37'/''T^/^6^F^' GeV. (173) 



The condition that the initial amplitude of the field ex- 
ceeds the critical value where the effects of gauge- and 
gravity-mediation on the effective potential are compa- 
rable is expressed as 



bo > W'^M^em^^^^GcV GeV. 



(174) 



Another condition for this situation to be realized is 
Mp ^ Teg, which can be rewritten as 



Mf > 5.0 X lO^f^Q^/^.-^/^^Tj 



l/4fl~l/16rTi5/4 
-RH,6 



GeV. 



(175) 



This implies that the allowed region in the parameter 
space {Q, Mf) can appear for low enough reheating tem- 
perature. It can be seen as foll ows. From the stability 
and B-DM conditions [Eqs.(|ll|) and ([l27|)], the largest 
possible Mp is 



Mf ~ 8.2 X 10^e^/^n]^^7js\ 



w'T^.pT GeV. 



(176) 



Therefore, comparing these two, we get the constraint on 
the reheating temperature as 

Trh < 1.5 X lO'£'/'ny'f3y'%R]{o'^;T G^V, (177) 

for generating the allowed region. On the other hand, 
we can also obtain the lower limit of Mf- In order to 
have successful nucleosynthesis, we need Trh ^10 MeV 
conservatively. We thus get the lower bound as Mf ^ 0.5 
GeV from Eq.( |l75|) . However, we will assume Mf ^100 
GeV, which may be the lowest possible scale for SUSY 
breaking, and this limit is severer. 

At the same time, we can estimate the largest possi- 
ble gravitino mass, taking into account the fact that the 



in general. 

The initial amplitude of the field is the last to be con- 
sidered. The initial value is determined by the balance 
of the Hubble mass and nonrenormalizable terms, so we 



have exactly same formulas as in Eq.(161). It depends 
only on the graviti no m ass as 00 ~ (7713/2^""'^)"'^^'"^^'' • 
Equating with Eq.(172), we get the relation between the 
gravitino mass and the reheating temperature. As easily 
seen, there is no allowed range for n = 4 and 5, while we 
have 



490 keV < TO3/2< 0.16 GeV, 
1.1 X 10^ GeV > Trh> 6.1 GeV, 



for n — 6, 



1.8 eV < m3/2< 0.16 GeV, 
1.1 X 10^ GeV > Trh> 72 MeV, 



for n = 7, 



6.4 X lO"^'"^ GeV < m3/2< 8.6 x 10"^ GeV 
1.1 X 10^ GeV > Trh> 10 MeV, 



(182) 



(183) 



(184) 



for n = 8, and so on. 

On the other hand, we can get the value of Mp required 
for both the Q-ball formation [Eq.(|l7l])] and the B-DM 
condition, in terms of the gravitino mass as 

Mf,, - 0.28ei/4f^;j/V7^/^yB.Yo'-;y'-3/2,GcV (185) 

Inserting the TO3/2 — Trh relation for each n into this 
equation, we have 

r 9.2 X 10-3f]i9/i2^7.f r^^/g^ GeV, (n = 6). 



Mf 



1.1 X iQ-^n 



-7o23/12^-l/6rp-5/3 



Q 



>,2 



?^Ki^,6 GeV, (n = 7). 



1.2 X lO-'^^%'m'^l^'T^l, GeV, 



(n = 



1.4 X lo-^^n 



17o31/12_-l/6.Ti-7/3 



Q 



m, 



TRti;GeY, (n = 9), 
(186) 



where we suppress the dependences on /3, 77s, and e, since 
they are trivial. These values have to satisfy the con- 
straint Eq.(175). They lead to the upper hmit for the re- 
heating temperature, and we obtain allowed regions only 
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f or n = 6, 7 . and 8, taking into account constraints (179), 
(180), and (181). Thus, we find the allowed region for the 
consistent scenarios as 



6.1 GeV < Tj^H< 29 GcV, 
8.3 X 10* GeV > Mp > 1.0 x 10"* GeV, 

0.16 GcV > m3/2> 7.0 x 10"^ GeV, (187) 

for n = 6, 

72 MeV < Trh< 1.9 GeV, 
8.3 X 10* GeV > Mp > 3.6 x 10^ GeV, 

0.16 GcV > m3/2> 4.5 x 10"^ GcV, (188) 

for n = 7, and 

10 MeV < Trh< 110 MeV, 
1.2 X 10* GeV > Mf> 10^ GeV, 

8.6 MeV > m3/2> 6.4 keV, (189) 

for n = 8, when we take — 100 GeV. 



X. DETECTION OF THE DARK MATTER 
Q-BALL 

A. Q balls in the specific logarithmic potential 

As is weh known, Q balls are stable against the de- 
cay into nucleons in the gauge mediation mechanism for 
SUSY breaking H,^. Therefore, they can be considered 
as the good candidate for the dark matter of the uni- 
verse. In addition, they can supply the baryon number 
of the universe. In the previous sections, we investigate if 
there is any consistent cosmological scenario for the dark 
matter Q ball with simultaneously supplying the baryon 
number of the universe. We have also seen the amount 
of the charge evaporated from a Q ball, so we can relate 
it to the amount of the dark matter. In this section, we 
see more conservative allowed region for the cosmologi- 
cal Q-ball scenario to work, and impose the experimental 
bounds in order to see if the scenario could exist. 

Provided that the initial charge of the Q ball is larger 
than the evaporated charge, we regard that the Q ball 
survives from evaporation, and contributes to the dark 
matter of the universe. This is expressed in Eq.(p3) 



Q < lO^^e 



23^3/2 / 



^ 7.4 



VTeV/ 



-12/11 



(190) 



Now we can relate the baryon number and the amount 
of the dark matter in the universe. As mentioned above, 
the baryon number of the universe should be explained 
by the amount of the charge evaporated from the Q balls, 
AQ, and the survived Q balls become the dark matter. 
This condition is Eq.(|7e^ 



VTcV/ ' 
IQ-i" and Qq < 1. 



(191) 



where we take rjB 

In order for the Q ball to be stable against the decay 
into nucleons, i.e., Eq/Q < 1 GeV, the following should 
be satisfied: 

Q>10^^(^)'. (192) 

For the usual gauge-mediation type of the Q ball, we 
obtain the allowed region for explaining the baryon num- 
ber of the unive rse, using three constraints, i.e., Eqs. 



(|190|),(|191D, and (|192| ). Figure |12| shows the allowed re- 
gion on {Q, 171(1,) plane with e = 1 in Eq.(191) ||l^. The 
shaded regions represent that this type of stable Q balls 
are created, and the baryon number of the universe can 
be explained by the mechanism mentioned above. Fur- 
thermore, this type of stable Q balls contribute crucially 
to the dark matter of the universe at present, if the Q 
balls have the charge given by the thick dashed line in the 
figure. Notice that this line denotes e = 1 case, and the 
allowed region in the parameter space will be narrower 
as e becomes smaller, and disappear for e ^ 10~^. 

As can be seen in Fig. |l^, several experi- 
ments constrain the parameter space. The Q ball 
can be detected through so-called Kusenko-Kuzmin- 
Shaposhnikov-Tinyakov (KKST) process. When nucle- 
ons collide with a Q ball, they enter the surface layer of 
the Q ball, and dissociate into quarks, which are con- 
verted into squarks. In this process, Q balls release ^ 1 
GeV energy per collision by emitting soft pions. This 
is the basis for the Q-ball detections H^,^. Lower left 
regions are excluded by the various experiments. The 
allowed charges are Q ~ 10^^ with = 100 GeV — 1 
TeV for e = 1, and future experiments such as the Tele- 
scope Array Project or the OWL-AIRWATCH detector 
may detect the dark matter Q balls. 

We also put the cons istent scenario in Fig. |l^, which we 
considered in Sec. |VIl| . Two thick lines inside the allowed 
region (shaded region) represent for the n = 5 and n = 6 
cases. Hatched line connected to the thick line for n = 5 
shows the allowed region if we do not assume that the A 
terms come from the nonrenormalizable superpotential. 
The current experiments exclude the n = 5 case, but the 
top edge of n = 6 case is allowed, which we show by 
circle in the figure. Although the n — 6 case sits at very 
interesting region, this dark matter Q ball may not be 
detected by future experiments as mentioned above. 

If we do not impose the condition that the evaporated 
charges from Q balls account for the baryons in the uni- 
verse, the only constraint is that the energy density of Q 
balls must not exceed the critical density. As mentioned 
in the previous section, this condition is Eq.(|7q): 



Q 9.3 xiO s KpSGevi 



/3 



6 X 10- 



8/5 



VTeV/ 



-12/5 
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FIG. 12. Summary of constraints on (Q, m^,) plain for the 
usual type of the Q ball. We also show the regions currently 
excluded by BAKSAN (B), Gyrlyand (G), and Kamiokande 
(K-1, K-2, K-3), and to be searched by the Telescope Ar- 
ray Project (TA) and OWL-AIRWATCH_(pA) in the future. 
For the details of experiments, see Ref. |^ . Two thick lines 
within the shaded region denote the consistent scenarios for 
n = 5 and n = 6 cases, respectively, and the only allowed con- 
sistent scenario is found inside the circle. Five thick dotted 
lines denote the constraints of the dark matter Q ball which 
do not account for the baryogenesis for Trh = 10^, 10"', 10*, 
10^ and 10'' GeV from the top to the bottom. 



The Q-ball charge depends on the reheating temperature, 
and thick dotted parallel lines show the constraints for 
Trh = 10^, 10^, 10'', 10^ and 10^ GeV from the top to 
the bottom in the figure. We can thus discover the dark 
matter Q balls in the future experiments, if the reheat- 
ing temperature is low enough. Notice that the lower 
reheating temperature is favored by the supersymmet- 
ric models because of evading the cosmological gravitino 
problem [[l9[. 



B. Q balls in the generic logarithmic potential 

Now we show the allowed region for the Q-ball forma- 
tion in the generic potential, but actual formation takes 
place when the thermal logarithmic term still dominates 
the potential. We show for the n = 6 case, which is the 
only successful scenario in this situation, in Fig. |l3|. The 
thick solid line denotes the natural consistent scenario, 
while the thick dashed line represents the allowed region 
on the B-DM condition line in a more general situation 
that the Q-ball formation occurs through the different 
mechanism. On these lines, the Q balls can account for 
both the dark matter and the baryons in the universe 
simultaneously. Most of the parameter space is excluded 
by currents experiments, but very interesting region such 
as Q ~ 10^'' and Mp ~ 5 x 10^ GeV is allowed. Notice 
that the dark matter Q balls with larger charges can be 
detected by the future experiments such as TA and OA, 
although they cannot play the role for the baryogenesis. 



C. Q balls in the gravity- mediation dominated 
potential 

Next, we move on to the 'new' type of Q ball. The 
constraints on parameter space (Q, '713/2) were obtained 
in Ref. Q. Here we only draw the results. The regions 
where both the amount of the baryon and the dark mat- 
ter can be explained are shown as thick lines in Fig. [T^ . 
Lower left regions are excluded by the various experi- 
ments. Notice that future experiments such as the Tele- 
scope Array Project or the OWL-AIRWATCH detector 
may detect this type of the dark matter Q balls, supply- 
ing the baryons in the universe, with an interesting grav- 
itino mass ^ 100 keV, if the the origin of the A terms 
differs from the nonrenormalizable superpotential. 

Although any consistent cosmological Q-ball scenario 
does not exist in the framework that both the initial am- 
plitude and the A terms of the AD field are determined 
by the nonrenormalizable superpotential, it is not nec- 
essarily true that the dark matter Q balls do not exist 
whether they can be the source for the baryons in the 
universe or not. In general, we have the constraint for 
the Q balls to be a crucial component of the dark matter. 
It reads as 



g ~ 3.5 X IQ^^VlQ 



6 X 10-^ 



VMeV/ 



VlQSGeVy 



(194) 



where Eq.(113) is used. 

Three dot-dashed lines are shown for Trh — 10'^, 10*^, 
and 10'', from the top to the bottom, respectively. There- 
fore, the dark matter Q balls can be detected in the future 
experiments for the low reheating temperature universe, 
which is free from the cosmological gravitino problem. 
Notice that the baryons have to be supplied by another 
mechanism in this case. 
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FIG. 13. Summary of constraints on the parameter space 
{Q,Mf) with n — 6 and — 100 GeV in the generic log- 
arithmic potential where the thermal terms are dominated 
when the Q-ball formation occurs. Dashed and dot-dashed 
lines represent the same current and future experiments, re- 
spectively, as shown in Fig. Lines (a), (b), and (c) denote 
the B-DM, survival, and stability conditions, respectively. 
Line (d) shows the conditions T > Mp. Notice that this 
condition depends on n. The thick solid line is the allowed 
region of the successful scenario for the dark matter Q balls 
supplying the baryons in the universe, while the thick dashed 
line represents more general cases. 
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FIG. 14. Restrictions by several experiments on (Q, 
plain for \K\ =0.01. We show the regions currently excluded 
by BAKSAN (B), Gyrlyand (G), and Kamiokande (K-1, K-2, 
K-3), and to be searched by the Telescope Array Project (TA) 
and OWL-AIRWATCH (OA) in the future. The thick lines 
represent for the gravity-mediation type of the Q ball to be 
both the dark matter and the source for the baryons of the 
universe for — 300 GeV, 1 TeV, and 3 TeV from the top 
to the bottom. Three dot-dashed lines denote the constraints 
of the dark matter Q ball which do not account for the baryo- 
genesis for Trh = 10^, 10^ and 10^ GeV from the top to the 
bottom. 



D. Delayed Q balls in the thermal logarithmic 
potential 

As shown in the previous sections, we can restrict the 
parameter space {Q,Mp) by several conditions for the 
Q-ball formation in the natural scenarios. In general, we 
have three conditions: the B-DM, survival, and stability 
conditions. Here, we write them down again: 



(a) Q ^ 3.2 X 10l^£3/2^3/2^^3/2^^-1^^^2 

(b) Q>1.2xl0i^ 



6' 



-8/ll,,r-4/ll 

',2 -'^-'f,6 ' 



(C) Q>10^H'I^ 



(195) 
(196) 
(197) 



We plot these lines in Figs. 15 and for n = 6 and 
n = 7 cases, respectively. In these figures, (a), (b), and 
(c) denote the B-DM, survival, and stability conditions. 



respectively, for fig = 1 and = 100 GeV. We also 
show the data from several experiments. These are the 
same as we have used above. The experimental lines are 
the same as that for the specific logarithmic potential, if 
we just replace by Mp. We plot these data, using the 
fact that the flux is related to the Q-ball mass, and the 
cross section is related to the Q-ball size. Therefore, these 
Q-ball parameters can be expressed in the same form as 
in the specific potential for the replacement of by 
Mp in the generic potential. Notice that the presence 
of the dark matter Q balls for — 1 TeV is almost 
excluded by experiments: Only Q balls with Q ~ 10^^ 
and Alp ~ 5.6 x 10"* GeV are allowed. 

Now we put the conditions for the Q-ball formation 
to take place in the natural scenarios. In Fig. [ij, ver- 
tical line (d) denotes for the condition that the thermal 
logarithmic potential dominates over the generic ones. 
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FIG. 15. Summary of constraints on the parameter space 
(Q, Mf) for the delayed- formed Q balls with the thermal log- 
arithmic term domination for = 100 GeV and n = 6. 
Dashed and dot-dashed lines represent the same current and 
future experiments, respectively, as shown in Fig. |l2| Lines 
(a), (b), and (c) denote the B-DM, survival, and stability 
conditions, respectively. Lines (d) and (e) show the condi- 
tions Teq ^ Mf and 7713/2 i$ 1 GeV, respectively. The thick 
solid line is the allowed region of the successful scenario for 
the dark matter Q balls supplying the baryons in the uni- 
verse. Thick dotted lines denote the formed Q-ball charges 
with Trh ~ 10^, 10^, and 29 GeV from the top to the bot- 
tom. 
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FIG. 16. Summary of constraints on the parameter space 
(Q, Mf) for the delayed-formed Q balls with the thermal log- 
arithmic term domination for = 100 GeV and n = 7. 
Dashed and dot-dashed lines represent the same current and 
future experiments, respectively, as shown in Fig. |l2| Lines 
(a), (b), and (c) denote the B-DM, survival, and stability con- 
ditions, respectively. Line (d) shows the condition Trh ^ 10 
MeV. The thick solid line is the allowed region of the success- 
ful scenario for the dark matter Q balls supplying the baryons 
in the universe. Thick dotted lines denote the formed Q-ball 
charges with Trh ~ 62 GeV and 10 MeV for the upper and 
lower lines, respectively. 



i.e., Teq ^ Mp. Another vertical line (e) represents for 
™3/2 ^ 1 GeV in the gauge- mediated SUSY breaking 
model. Three thick dotted lines show the charges of the 
produced Q balls in the natural scenarios with the re- 
heating temperature 1.0 x 10^, 1.0 x 10"^, and 29 GeV 
from the top to the bottom, respectively. The thick solid 
line is the final answer, which represents the dark matter 
Q balls supplying the baryons in the universe in the natu- 
ral scenario. Two remarks are following. (1) Those dark 
matter Q balls which do not account for the baryogenesis 
may be detected by the Telescope Array Project (TA), if 
their charges are ~ 10^^ and the reheating temperature is 
^ 10^ GeV. (2) There is no cosmological gravitino prob- 
lem for this scenario. 



In Fig. 16, line (d) denotes the condition that the re- 
heating temperature should be higher than 10 MeV in 
order for the nucleosynthesis to take place successfully. 
We also assume Mp ^ 100 GeV. Two thick dotted lines 



represent the Q-ball charges from the formation mech- 
anism with the reheating temperature 62 GeV and 10 
MeV for upper and lower lines, respectively. Thick solid 
line shows the successful region that the Q balls account 
for both the dark matter and the source for the baryons 
in the universe simultaneously. Most of this region is al- 
lowed experimentally, and, in particular, the dark matter 
Q balls with Q ~ 10^^ may be detected by the TA in the 
future. Notice again that cosmological gravitino problem 
is avoided for such low reheating temperature. 

Now we will mention the specific logarithmic poten- 
tial, where Mp = rn^. We impose three general con- 
ditions; namely, the B-DM, survival, and stability con- 
ditions. They are written as [cf. Eqs.(191), ( |190| ), and 

(198) 
(199) 



(192)] 
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FIG. 17. Summary of constraints on the parameter space 
{Q,m^) for the delayed-formed Q balls with the thermal log- 
arithmic term domination for n = Q. Experimental lines are 
the same. Lines (a), (b), and (c) denote the B-DM, survival, 
and stability conditions, respectively. Thick dotted lines de- 
note the formed Q-ball charges for Trh — 10^, 3 x 10*, and 
5 X 10^ GeV from the top to the bottom. The thick dashed line 
represents the condition Teq ^ m^. The thick solid line is the 
allowed region of the successful scenario for the dark matter 
Q balls supplying the baryons in the universe. Solid vertical 
lines denote the gravitino mass ranges, 7713/2 = 1 GeV, 100 
MeV, 10 MeV, and 1 MeV from the left to the right. 
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FIG. 18. Summary of constraints on the parameter space 
(Q,7n0) for the delayed-formed Q balls with the thermal log- 
arithmic term domination for n = 7. Experimental lines are 
the same. Lines (a), (b), and (c) denote the B-DM, survival, 
and stability conditions, respectively. Thick dotted lines de- 
note the formed Q-ball charges for Trh — 10'' and 80 GeV for 
the upper and lower lines. The thick dashed line represents 
the condition Teq ^ r7i^. The thick solid line is the allowed 
region of the successful scenario for the dark matter Q balls 
supplying the baryons in the universe. Solid vertical lines de- 
note the gravitino mass ranges, 7773/2 = 100 keV (left) and 10 
keV (right). 



(c) Q>10'^m' 



sTcV 



(200) 



which are shown in Fig s. and |l|. 

The charge of the formed Q-ball and the initial ampli- 
tude of the field are expressed as 



1)0 5.3 X wriQ^'y^ 



26n4/5fl4/5_4/5 ^4 -4 (r,r..-. 

.14o-2/5fl-l/10 



g ^ 1.7 X 10"^r!Q'"/3;'^m:';;„,,T^„ ,m 



^ "^/toV^Ih, 5 '^3/2 , Gov ' 



(202) 



respectively. The initial amplitude is determined by the 
balance between the Hubble mass and nonrenormalizable 



terms, and combining it with Eq.(202), we get the rela- 
tion between tti^ and 7773/2 as 

771^ ~ 4.9 X lO^riQ/^y^rngT^^Lv GeV, (203) 

for 77 = 6. We thus get the charge of the produced Q 
balls: 



Q - 1.9 X 



Q 



RH,5' 



36/5 
^(A,TcV' 



(204) 



which we plot for Trh = 10^ 3 x 10'', and 5 x 10^ 
GeV from the top to the bottom by thick dotted lines 
in Fig. ^ In order for this situation to take place, 
Teq 7770 , which Icads to the constraint for the reheating 
temperature in terms of 777^ as 



Trh > 4.4f7i/'/3^ 



'^''^'/'"-'/tcV GeV. 



Inserting this constraint into Eq.(201), we have 



Q > 6.8 X 10^°n, 



-3/5 52/5 



.,ToV' 



(205) 



(206) 



which is shown by the thick dashed line in Fig. As 
is seen in the figure, the experimentally allowed region 
is very narrow, and the consistent natural scenario will 



work only for Q ^ 10"''', 777^ ~ 5 x 10 GeV, 7773/2 ^ 1 
GeV, and Trh ~ 3 x lO"* GeV. Notice that there is no 
cosmological gravitino problem also in this case. 
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On the other hand, the relation be twee n and TO3/2, 
which is obtained by equating Eg. ( ^02| ) and the value 
determined by the balance between the Hubble mass and 
the nonrcnormalizable term, is written as 



35 



GeV, 



(207) 



for n — 7. We thus get the charge of the produced Q 
balls: 



Q - 1.9 X 



28o-28/5fl-4/5^4 



A" 



RH.5' 



36/5 
^<*,TcV' 



(208) 



where we plot for Tj^h ~ 10'' and 80 GeV for the upper 
and the lower thick dotted lines, respectively, in Fig. 1^ . 
Taking into account the condition 



eq 



> 



we obtain 



(209) 



where Eq.( |205| ) is used. We plot this constraint by the 
thick dashed line in Fig. The allowed region for the 
dark matter Q balls, which are also the source for the 
baryons in the universe, is shown by thick solid line. This 
region does not suffer from the current experimental limit 
and this dark matter Q ball with Q ~ 10^^ and TO3/2 ^ 
100 keV may be detected by TA in the future. Notice 
that the reheating temperature is low enough to avoid the 
cosmological gravitino problem. If we abandon to explain 
the baryons in the universe by the charge evaporation 
from Q balls, the dark matter Q balls may be found by 
TA and also OWL-AIRWATCH in the future. 



E. Delayed Q balls in the generic logarithmic 
potential 

As mentioned in the previous section, we have obtained 
the consistent scenarios of the dark matter Q ball explain- 
ing the baryons of the universe naturally, with the initial 
conditions determined by the balance of the Hubble and 
the nonrenormalizable terms in the effective potential, 
with n = 6, 7, and 8. However, the parameter space 
is restricted by the current experiments. We plot the 
allowed regions and experimentally excluded regions for 
n = 6, 7, and 8 in Figs. and |2^, respectively. 

In these figures, lines (a), (b), and (c) are the B-DM, 
survival, and stability conditions, respectively. Lines (d) 
m Figs. |l| and |^ are the condition Mp ^^J^eg, while 
Tn,H ^ 10 MeV is shown as Hne (d) in Fig. Lines (e) 
in Figs. |l9| and ^ are just the upper limit for Mp. 

We can see some parameter regions which are experi- 
mentally allowed in each figures. Typically, there are two 
types: One is for Q - lO^^ and Mp - 5 x 10'* GeV, and 
the other is for Q - lO^^ and Mf - 3 x 10^ GeV. In the 
latter case, the dark matter Q balls may be detected in 
the future experiments. Notice that the required reheat- 
ing temperature is rather too low (100 MeV — 10 GeV), 
which is very difhcult to realize in the actual inflation 
model. 
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FIG. 19. Summary of constraints on the parameter space 
(Q, AIf) for the delayed- formed Q balls with the generic log- 
arithmic term domination for n — 6 and = 100 GeV. 
Experimental lines are the same. Lines (a), (b), and (c) 
denote the B-DM, survival, and stability conditions, respec- 
tively. Line (d) represents the condition Mp ^ Teq, and line 
(e) is just the upper limit for Mp. 



XI. CONCLUSION 

We have investigated thoroughly the Q-ball cosmology 
(the baryogenesis and the dark matter) in the gauge- 
mediated SUSY breaking model. Taking into account 
thermal effects, the shape of the effective potential has 
to be altered somewhat, but most of the features of the Q- 
ball formation derived at zero temperature can be applied 
to the finite temperature case with appropriate rescal- 
ings. We have thus found that Q balls are actually formed 
through Affleck-Dine mechanism in the early universe. 

We have sought for the consistent scenario for the dark 
matter Q ball, which also provides the baryon number of 
the universe simultaneously. For the consistent scenario, 
we adopt the nonrcnormalizable superpotential in order 
to naturally give the initial amplitude of the AD field and 
the source for the field rotation due to the A-term. As op- 
posed to our expectation, very narrow parameter region 
could be useful for the scenario in the situations that the 
Q balls arc produced just after the baryon number cre- 
ation. In addition, we have seen that current experiments 
have already excluded most of the successful parameter 
regions. 

Of course, if the A-terms and/or the initial amplitude 
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FIG. 20. Summary of constraints on the parameter space 
{Q, Mf) for the delayed- formed Q baUs with the generic log- 
arithmic term domination for n — 7 and — 100 GeV. 
Experimental lines are the same. Lines (a), (b), and (c) 
denote the B-DM, survival, and stability conditions, respec- 
tively. Line (d) represents the condition Mf ^ Teq, and line 
(e) is just the upper limit for Mf- 
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FIG. 21. Summary of constraints on the parameter space 
(Q, Mf) for the delayed- formed Q balls with the generic log- 
arithmic term domination for n = 8 and = 100 GeV. 
Experimental lines are the same. Lines (a), (b), and (c) 
denote the B-DM, survival, and stability conditions, respec- 
tively. Line (d) represents the condition Trh ^ 10 MeV. 



of the AD field are determined by other mechanism, the 
cosmological Q-ball scenario may work. Then, the stable 
dark matter Q balls supplying the baryons play a crucial 
role in the universe. 

We have also found the new situations that the Q-ball 
formation takes place when the amplitude of the fields 
becomes small enough to be in the logarithmic terms in 
the potential, while the fields starts its rotation at larger 
amplitudes where the effective potential is dominated by 
the gravity-mediation term with 'positive K-term. This 
allows to produce Q balls with smaller charges while cre- 
ating larger baryon numbers. In this situation, there is 
wider allowed regions for naturally consistent scenario, 
although the current experiments exclude most of the 
parameter space. Notice also that rather too low reheat- 
ing temperature is necessary for larger n scenario to work 
naturally. This aspect is good for evading the cosmolog- 
ical gravitino problem, while it is difficult to construct 
the actual inflation mechanism to get such low reheating 
temperatures, in spite of the fact that the nucleosynthesis 
can take place successfully for the reheating temperature 
higher than at least 10 MeV. 

So far we have investigated the possibility that the dark 
matter Q balls which supply the baryon number of the 



universe could be really produced in the naturally con- 
sistent scenario, where the initial amplitude of the AD 
field and the A terms for the field rotation are deter- 
mined by the nonrenormalizable superpotential. Here, 
we have explicitly imposed that the produced Q balls 
must be stable against the decay into nucleons in order 
for Q balls themselves to be the dark matter of the uni- 
verse. However, there may be such situations that Q balls 
decay into lighter particles when the stability condition 
is broken. This situation is almost the same as the usual 
AD baryogenesis when we consider only the baryons in 
the universe. Moreover, if the decay products include 
the lightest supersymmetric particles (LSPs), they may 
account for the dark matter of the universe. |^ 

If this is the case, gravitino may be the LSP to be the 
dark matter. The number density of LSP depends on 
the lifetime of the Q ball, as investigated in Ref. ||l^ for 
the gravity-mediated SUSY breaking model. If it decays 



^ We thank M. Fujii for indicating the situation that the Q 
balls can decay and become the source for the baryons and 
the LSP dark matter simultaneously in the gauge-mediated 
SUSY breaking model. 
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before the temperature drops the freeze-out temperature 
of LSPs, the density follows the thermal value. On the 
other hand, when its lifetime is long enough, the ratio 
of the number of the baryons to that of LSPs is fixed, 
and we can directly estimate the relation between the 
amounts of the baryons and the LSP dark matter of the 
universe. 
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